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Abstract 

Recently, the development of chiral perturbation theory has allowed 
the generation of rigorous low-energy theorems for various hadronic pro¬ 
cesses based only on the chiral invariance of the underlying QCD La- 
grangian. Such techniques are highly developed and well-tested in the 
domain of pionic and kaonic reactions. In this note we point out that 
with the addition of a few additional and reasonable assumptions similar 
predictive power is available for processes involving the eta meson. 
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1 Introduction 


It has long been the holy grail for particle and nuclear knights to generate 
rigorous predictions from the Lagrangian of QCD 

>Cqcd = - - m)q (1) 

where 

G — 3^-A.y 3yA^ 

D^q = {df, - igA^)q. (2) 

However, despite the ease with which one can write this equation, because of its 
inherent nonlinearity and the large value of the coupling constant—5^/47r ~ 1— 
progress in this regard has been slow. One approach—lattice gauge theory— 
holds great promise, but is currently limited by the need for large computational 
facilities 1^. A second tack, that of perturbative QCD, exploits the feature of 
asymptotic freedom—the vanishing of the running coupling constant at high 
momentum transfer]^. However, such predictions are valid only for the very 
highest energy processes. It is gratifying then to see that in recent years a third 
procedure has become available, that of chiral perturbation theory (yPT) which 
exploits the chiral symmetry of QCD and allows rigorous predictive power in 
the case of low energy reactions. This technique, based on a suggestion due 
to Weinberg^, was developed (at one loop level) during the last decade in 
an important series of papers by Gasser and Leutwyler and others Q. The 
idea is based on the feature that the QCD Lagrangian—Eq. |I] —possesses a 
global SU{3)l X SU{3)ji (chiral) invariance in the limit of vanishing quark 
mass. Such invariance is manifested in the real world not in the conventional 
Wigner-Weyl fashion. Rather, it is spontaneously broken, resulting in eight 
light pseudoscalar Nambu-Goldstone bosons— tt, K, r] —which would be massless 
if the corresponding quark masses also vanished |^. While the identihcation of 
this symmetry is apparent in terms of quark/gluon degrees of freedom, it is not 
so simple to understand the implications of chiral invariance in the arena of 
experimental meson/baryon interactions. 

Early attempts in this direction were based on current algebra/PCAC 
methodsl^, yielding relationships between processes differing in the number 
of pions, e.g. 

lim<H7r“|0|A>=g <H|[E“,0]|A> (3) 

where = 92.4 MeV is the pion decay constant 0. However, it was soon 
realized that the most succinct way to present these restrictions is in terms of 
an effective chiral Lagrangian, the simplest (two-derivative) form of which is, in 
the Goldstone sector [||, 

4^) = ^TrDi^UDf^U^ + ^TT2Bom{U + U^) + --- (4) 
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where 


U = exp 



(5) 


is a nonlinear function of the pseudoscalar helds, m = 
quark mass matrix, 


2Bo = 


2m|. 
TO„ + rus 


2ml 

rriu + rud 


(m„,TOd,TOs)diag is the 


( 6 ) 


is a phenomenological constant, is the covariant 

derivative coupling to external fields and F is the pion decay constant 

in the limit of chiral symmetry. Although these are only two of an infinite 
number of terms, already at this level there exists predictive power-e.g., tree 
level evaluation of yields the familiar Weinberg predictions (at to2)) 

for S-wave tt — tt scattering lengths]^ 


a 


0 _ 
0 — 


7m-j^ 

3^ 


On = 


167rF2 


(7) 


which are approximately borne out experimentally. Loop diagrams, of course, 
produce terms of ,m'^) and contain divergences. However, just as in 

QED such infinities can be absorbed into renormalizing phenomenological chiral 
couplings, and the most general ’’four-derivative” Lagrangian has been given by 
Gasser and Leutwyler[^ 

= L,{TtD^UD^U^)^+L2{TyD^UD^U^)^ 

+L3TT{Df,UDf^U'<f + L4TrD^Z/L)'"C/t)Tr(m(C/ -k G^)) 
+L3TTDf,UD>^U'fmiU + C/t) L6(Trm(t7 -k 
-kL7(Trm(C/ - -k LsTrimUmU + 

+iLQ{TiF/;;^Df^UD''U'< + TtF^^D>^UD’'U'^) 

+LioTrF^^^C7F«'^"Z/t (g) 

where Fj;^, F^^ are external held strength tensors dehned via 


pL,R 



d^F^'^ - d,F^’^ - F^-^] 

Vf,+A^ F/^ = Vf,-Af,. 


(9) 


Here the bare Li coefficients are themselves unphysical and are related to em¬ 
pirical quantities L”(/i) measured at scale p, via 


T'i (a^) — Li + 


T, 


327r2 


1 477 

- -I- In^ -1-1-7 
e 


( 10 ) 


where are constants dehned in ref. 4b and e = 4 — d is the usual parameter 
arising in dimensional regularization, with d being the number of dimensions. 


2 







L\ 

T ^ 

rr 

-^5 

T ^ 

rr 

-^10 

0.71 ±0.28 

2.01 ± 0.37 

2.7 ± 0.3 

7.7 ± 0.2 

- 5.2 ± 0.3 


Table 1: Empirical values of Chiral Expansion Parameters(x 10 with /i = 


Gasser and Leutwyler have obtained empirical values for the phenomenological 
constants values for some of which are given in Table 1. 

At the four-derivative level it is also necessary to include the contribution 
of the anomaly, which in the case of coupling to the photon field has the 
form|p^ 


^wzw{U, Afj_) = rwzw(G) 

+ d^x[eA^Tr{Q{R,R^Rfs + L.L^L^)) 

- ie^Ff^^A^Ti {q^{Rp + L0) + ^{QU^QUR/s + QUQU^L^)^ 


( 11 ) 


where = [d^U'^)U,L^ = Udf^W and Twzw is independent of the photon 
field and will not be needed for our purposes. A corresponding form involving 
coupling to a general nonabelian gauge field can also be written, but is lengthy 
and will not be given here|^. 

In a series of recent papers it has been conclusively demonstrated that this 
chiral effective action formalism provides a succinct and successful description 
of low energy electroweak interactions of pions and kaons[^. Specifically, the 
reactions given in Table 2 are successfully described in terms of GL parameters 
Lio(/r). Glearly, there are far more reactions than parameters, and this 
overdetermination enables one to construct required relationships between ex¬ 
perimental quantities, the empirical validity of which constitutes a strong test of 
the chiral formalism and indeed thereby of QCD itself|]^ . Such tests are found 
to be well satisfied, with the possible exception of the relationship between the 
axial structure function in radiative pion beta decay— Ha —and the charged pion 
polarizability—However, recent work indicates that this may not be 
a problem and in any case a number of experimental efforts are underway to 
retest this critical strict ure[p^. 

In the pion and kaon arena then one finds strong evidence for the correctness 
and utility of chiral perturbative techniques, and it is an obvious next step to 
attempt to extend this success into the eta sector, which is the subject of this 
note. This examination of the eta system is important both as a theoretical 
exercise and because of the existence of high intensity sources of etas|^. In 
the next section then we study the ability of the chirally inspired methods 
to make reliable predictions for eta decay processes. We examine only the 
’’allowed” modes —77 —> 27 , Stt, 27r7, Tr2^, STry, i.e. those modes which can occur 
assuming only isospin violation or the anomaly, eschewing the temptation to 
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Reaction 

Quantity 

Theory 

Experiment 

1 

+ 

T 

< rj >(fm^) 

0.44“ 

0.44 ±0.02 

7 ^ K+K- 

A 

V 

0.44 

0.34 ±0.05 

TT"*' ^ 7r+i/e7 

V(to^^) 

0.027 

0.029 ±0.017 


hA/hv 

0.46“ 

0.46 ±0.08 


{hv + hA){m~^) 

0.038 

0.043 ±0.003 

TT'^ —s- 

TAlhy 

2.6 

2.3 ±0.6 

77r+ — ■jTr+ 

{as -l-/3M)(10“'‘(fm^) 

0 

1.4±3.1 


a£:(10“‘*(fm^) 

2.8 

6.8 ± 1.4 

77 ^ TTTT 

a£:(10“‘^(fm^) 

2.8 

2.2 ± 1.6 

K —>■ Tre'^iye 

e = /-(o)//+(o) 

-0.13 

-0.20 ±0.08 


A+(fm^) 

0.067 

0.065 ±0.005 


0 

to 

0.040 

0.050 ±0.012 


Table 2: Chiral predictions and data in the radiative complex of transitions. 
The superscript a indicates that this parameter was used as input and is not a 
predicted quantity. 


analyze “rare” processes such as rj ^ 87 , as these have been well-discussed 
elsewherefp^. Finally, we summarize our results in a concluding section III. 


2 Eta Decay Processes 

Strictly from a kinematic perspective, inclusion of the ? 7 ( 547 ) as part of the 
chiral formalism should not be a problem, as the eta and kaon are roughly 
degenerate in mass, and as mentioned above the kaon (and pion) predictions 
obtained in this way are quite successful. Rather the real challenge involves 
mixing with f]'{958), which lies outside the simple chiral SU{3)l x SU{3)ii 
framework. To lowest order things are simple—in the chiral limit the pseu¬ 
doscalar mass spectrum would consist of a massless octet of Goldstone bosons 
plus a massive SU{3) singlet {rjo). With the breaking of chiral invariance the 
octet pseudoscalar masses become nonzero and are related, at first order in 
chiral symmetry breaking, by the Gell-Mann-Okubo formula 

« (0.57 GeV)^ (12) 

where r/g is the eighth member of the octet. At this same order in symme¬ 
try breaking the singlet ryo will in general mix with ryg producing the physical 
eigenstates 7 , 77 ' given by 


77 = cos 0ris — sin Orjo 

■q' = sin 077g-|-cos 0770. (13) 
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The mixing angle 9 can be determined via diagonalization of the mass matrix 
(written in the 773 , ?7o basis) 


II 

IM 

to 

(M pro 

"l08 

0 


V "^08 

<0 


(14) 


Here rri^^ is given in Eq. while TOQgjTO^^, and 6 are unknown. Diagonalizing 
and fitting these parameters with the two known masses yields the results 


9 = -9.4°, m^8 = -0.44m|:, = 0.95 GeV (15) 


However, there is good reason not to trust this simple and lowest order analysis, 
since higher order chiral symmetry breaking terms can generate important mod¬ 
ifications. For example, taking the leading log correction arising from Figure 1, 
we find[p^ 


m 


2 

Vs 


4 2 1 2 

(0.61 GeV)2 


_ 2 m|- 

3 (4^F,)2 
for 77 = 1 GeV, 


(16) 


for which diagonalization of the mass matrix yields 

6l«-19.5°, = =-0.90 GeV (17) 

suggesting a doubling of the mixing angle. Of course, this is just an approximate 



Figure 1: Mass and wavefunction renormalization diagram. 

result. However, a full one loop calculation using yPT yields essentially the same 
result Q, and consequently in our analysis below we shall use the value in Eq. 
1^, i.e. 

1 2 /2 

sin0«—- cos0«——. (18) 

o o 

It is also intriguing that this solution is consistent with the assumptions of simple 
U(3) invariance wherein 778,770 have the same quark wavefunction, leading to 

TTlQg 2^/2 / 771 — 777s 

m\ 3 \iJi + ms 
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At this same (one-loop) level of symmetry breaking there is generated a shift 
in the lowest order value of the pseudoscalar decay constant Fp. Thus one finds 
from the diagrams in Figure 2, in leading log approximation, 




Fr,. 


F 

F 


^ '"-K 

2 (47rF^)2 

3 to|- 

2 (47rF^)2 


In 


In 


m 


K 


O 

m 


K 


1.12F 

1.25F^ for GeV. (20) 


Once again these estimates are in excellent agreement with those given in the 
full one-loop analysis Q. (We note in addition that the corresponding prediction 


Fk 1 ml ml 3 

F^ 4(47rF,)2 ^2 8 (47rA,)2 ^2 


1.22 


is quite consistent with the experimental value 1.22 ± 0 . 01 ) 


( 21 ) 




Figure 2: Loop diagrams leading to renormalization of the pseudoscalar decay 
constant. Here the symbol x indicates coupling to the axial current. 


2.1 r; —>■ 77 

With this introductory material in hand we can now confront the subject of our 
report—that of eta decay. First consider the dominant two-photon decay mode, 
which to leading order arises due to the anomaly. In the analogous 7 r° —> 77 
case we find from Eq. pd| 

Amp = 1 ^ 77 (0)e'"'"“^e^fc^e'„fc^ with E^^“(0) =-^^ = 0.025 GeV^ 

( 22 ) 

General theorems guarantee that this result is not altered in higher orders of 
chiral symmetry breaking]^ and, using the experimental value 

r(7r° ^ 77 ) = (7.7±0.7) eV, (23) 

we determine 

= (0.0250 ± 0.0005) GeV^ (24) 
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in excellent agreement with the theoretically predicted value and eloquently 
confirming the value Nc = 3 as the number of colors. Strictly speaking, the 
prediction of the anomalous four-derivative Lagrangian Eq. should be in 
terms of F rather than Indeed the difference between the value given in 
Eq. ^ and the strict four-derivative prediction involves terms of dimension six 
and is higher order in the chiral expansion. Clearly, however, our prediction for 
^.,^ 77 ( 0 ), which arises from what we shall term extended-yPT, is in excellent 
agreement with experiment. Nevertheless, although very reasonable, this is not 
a firm prediction of yPT itself. 

The ri,r]' ^ 77 couplings also arise from the anomalous component of the 
effective chiral Lagrangian and, in the extended yPT approximation, have the 
values 

F„77(0) = (^^cos9-2V2^sm0^ 

F^'77(0) = ^^ 77 ( 0 ) f ^ Q 2a/ 2^ cos 6*)^ . (25) 

v3 \i'8 ho / 

Using the experimental numbers 

r{rj 77 ) = (0.51 ± 0.05) keV T{r]' 77 ) = (4.7 ± 0.7) keV (26) 

we find 


F^.^^(0) = 0.0249 ± 0.0010 GeV^ h'r,'77(0) = 0.0328 ± 0.0024 GeV^ 

(27) 

In order to solve this system, we require an additional assumption since there are 
three unknowns— Fq, Eg, 6 —but only two pieces of data—Eq. |^. The standard 
approach at this point has been to use the leading log prediction from one-loop 
chiral perturbation theory— 


F^ 


= 1 - 


m 


K 


{AnF^y 


■In 




ml 


(47rF,)- 


■ In 


to: 




-1.30 at /r~lGeV (28) 


as input, and then to solve for Fq,9, yielding 

6»--20°, — w 1.04 


(29) 


It is intriguing that these results from two phton decay are quite compatable 
with those obtained from the one-loop analysis of the mass matrix— i.e., 0 — 
— 20 ° and Fq/Ft^ consistent with the value of unity which one would obtain if 
the singlet state and the pion were to have the same quark wavefunction. 

Closely related to the above modes are the associated Dalitz decays— 
n,r],r]' ye+e”—which have been studied at DESY|^. These are tradi¬ 
tionally parameterized in terms of a dipole slope parameter Ap such that 

^^ = (1 + ^)^ with s={p++p-f. (30) 


7 






The experimentally obtained values 


= 0.75 ± 0.03 GeV = 0.84 ± 0.06 GeV A^, = 0.79 ± 0.04 GeV 

(31) 

are in reasonable agreement with the vector dominance predictions [Em 


K = 


A^ = m 


A^, = m 


« (0.77 GeV)^ 

2 / 3cos0 — 6-\/2sin0 

p.“ ^5 _ 2^2^ COS0 — (5 + ^^)-\/2sin0^ 

^ 3 sin 0 + 6-^/2 cos 9 


(5 — 2^2) sin0 + (5 + ^‘^)\/2cos9^ 
where = {ml /ml) 


(0.75 GeV)2 
(0.83 GeV)2 


(32) 


These results are consistent with the observation that vector/axial dominance 
provides a remarkably successful representation of the values of the GL param¬ 
eters Ll{^) obtained empirically]^ and suggest an additional extension of our 
(already) extended yPT to include an effective vector-dominated Lagrangianp^ 


with 


and 


>Ceff 

= £vpp + £^7 + £vvp + Cppp^ 

(33) 

£vpp 

= ^TrH^[</>,9'^</>] 



~ -I- 


£vvp 

= -^5yvpe'‘"“^Tr(a^l49c.H/3</>) 


£pPP7 


(34) 


35^ 

(35) 


in order to understand the six-derivative component of the chiral expansion. 
Here g is given by the KSRF relation as® 


y 2p2 


In this formalism then the amplitude for tt^ 
diagram in Figure 3 to be 


2e 


(36) 

77 is determined from the 


F^77(0) - 3 2 2 (2ep;ff) ff-P- - 


(37) 
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as the value required by the anomaly. 

With this background in hand, we can now discuss a second important de¬ 
cay mode of the eta—that of 77 ^ TTTry which also arises from the anomalous 
component of the Lagrangian. 
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Figure 3: Vector dominance diagram responsible for —> 77 . 


2.2 T] —>■ 7r~^7r 7 

In the previous section we performed an analysis of the QCD anomaly as man¬ 
ifested in the two photon decay of the pseudoscalar mesons and, by use of the 
one loop leading log value for Fg/F^, were able to determine a solution for the 
77,77' mixing angle which is close to that found in the mass analysis together 
with a value for Fq/F^. which is near that which results from the assumption 
that the 770 and pion have the same wavefunction. While this is somewhat satis¬ 
fying, it is intriguing to inquire whether one can assess the mixing angle purely 
phenomenologically. We shall show below that this question can be answered in 
the affirmative, provided one utilizes the additional information available in the 
anomalous decays 77,77' —> 7r+7r“7. Such processes involving a photon coupled 
to three pseudoscalar mesons involve the anomaly and, at zero four-momentum, 
are completely determined from Eq. |^. However, inclusion of higher order ef¬ 
fects generates structure and study of such processes requires proper attention 
to issues of unitarity and final state interactions. Before considering 77,77' decay, 
however, we consider first the closely related case of 7 ^ 7r+7r“7r°. At zero 
four-momentum the anomaly requires 

Amp(37r-7) = A(s+_, s+o, s_o)e''''“^e;,p+,.73-aPo/3 

eN 

where A( 0 , 0 , 0 ) = i 2 ^ 2 p 3 ^ GeV~^ and Stj = {pt + pj)'^ (38) 

and one might suspect that vector dominance could reproduce this result di¬ 
rectly, as in the case of 7r° ^ 77 discussed above. However, this turns out not 
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to be the case. Rather, use of the diagram shown in Figure 4a yields 


^(s, t, u) 


eiVc f ml ml 

2A'k'^F^ yml — s ml — t 



(39) 


which at zero four-momentum is 50% larger than the value given by the anomaly. 
The resolution of this problem is well-known and arises from the a direct 7 — Stt 
coupling—Figure 4b—given in £ppp.^ whose origin presumably is from unspec¬ 
ified high-momentum-scale processes . Addition of this contribution to the 
7 — Stt process yields an amplitude 


A(s, t, u) 


eNc I ^ ^ ^ ^ ^ 

127r^F’3 2 \ml — s ml — t ml — u) 


(40) 


which has the structure required by vector dominance, but also agrees with the 
value required by the chiral anomaly at zero four-momentum. Inclusion of loop 
corrections modifies Eq. ^ slightly but does not change the general form. 




Figure 4: Vector dominance diagrams responsible for the anomalous process 

77r ^ TTTT. 


The 7—Stt reaction has been studied experimentally via pion pair production 
by the pion in the nuclear Coulomb field and yields a number 

A(0,0,0)exp = 12.9 ± 0.9 ± 0.5 GeV^ (41) 

in apparent disagreement with Eq. ^ and suggesting the value Nc « 4! This 
value was obtained, however, assuming no energy dependence of the amplitude 
and is reduced to 


A(0,0,0)exp = 11.9 ± 0.9 ± 0.5 GeV^ (42) 

if Eq. 1 is utilized, but is still too large. The most likely conclusion is that 
this an experimental problem associated with this difhcult-to-measure process, 
but it has recently been pointed out by Ametller, Knecht, and Talavera that 
an important electromagnetic effect—the photon exchange diagram connecting 
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7r*^77* s-nd 7 r+ 7 r “7 vertices—can reduce this number by another 1 x 10 “^ or 
so[^. In any case a new high-precision experiment is clearly called for, and 
this has been accomplished at JLab using the CLAS detector and the reaction 
7 P ^ 7 r+ 7 r“n. Such a measurement has been also been proposed at Da<i>nep^. 

The JLab experiment is currently being analyzed and this high statistics 
measurement will require an equally careful theoretical analysis in order to pro¬ 
duce the desired extrapolation from the rho resonance region, where most of the 
data has been obtained, to the zero four-momentum point where the anomaly 
stricture obtains. This issue has been approached in a number of authors: 

i) Holstein has used a simple matching of the one loop chiral correction to 
the rho dominance formj^. 

ii) Truong has utilized a unitarization based upon the Omnes-Muskhelishvili 
method Q. 

iii) Hannah unitarizes the amplitude using the inverse amplitude 
procedure]^. 

Regardless of the method used, the results are similar and are somewhat robust. 
The resulting value of the anomaly obtained in a preliminary analysis of the 
CLAS measurement are consistent, within a signihcant uncertainty, with the 
expected number—three. However, a definitive value awaits further analysis. 

Having warmed up on the 7 — Stt process, it is now straightforward to con¬ 
struct the analogous rj —>■ 7 r+ 7 r “7 amplitude. Using the extended xPT assump¬ 
tion we findj^ 

Amp( 7 ? ^ 7 r+ 7 r“ 7 ) = B{s+-,s+j, S-^)e^''°'^e*^p+^p-ak^p (43) 

with the anomaly stricture yielding 

However, the physical region for the decay—4m^ < St^tv < —is far from 
the zero-momentum point which is constrained by the anomaly. One indi¬ 
cation of this fact is that the decay rate obtained via neglect of momentum 
dependence—r^'^l* 7 r 7 r 7 = 35.7 eV—is signihcantly different from the experimental 
value—= 64± 6 eV. For the r]' channel the situation is, of course, much 
worse. The experimental value—= 61 ± 5 keV—is a factor of twenty 

larger than the value F*'®^ , _ =3 KeV obtained via use of the simple anomaly 

prediction 

(0,0,0) =-- f — sin0-I-V^—cosfl') (45) 

Thus proper inclusion of momentum dependence is essential. The ij —>■ 7 r^ 7 r “7 
spectrum was measured in the experiment of Gormley et al.|4l|] and was found 
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to be approximately fit in terms of a pure (width-modified) p-dominated ma¬ 
trix element. This result is not in agreement, however, with the simple vector 
dominance prediction— cf. Eq. ^0| —which would require 

^^)ltheo ~ 1 + 3—J + (46) 

and corresponds instead to 

1 + 2^+ ■■■. (47) 

nip 


Thus a careful look at the unitarization procedure is called for. 

We begin by noting that a one-loop chiral perturbation theory calculation 
gives 


Bli ‘°°P(s,s^,r) = ■B^(0,0)[1-f ^^-^((-4m2-h is^^)ln^ 


+ 


20 

y 


2m2‘ 


(48) 


where 


F{s) = 




-H- 


— 2 s > Ami 


2(1-4;^) 


4m2 —s . _1 

—^— tan 


Am'^—s 


— 2 s < Am 


— TT 


(49) 


while the vector dominance picture {cf. Figure 5) yields 


B^{s^^)=B^{0,0,0)(l + ^ ) (50) 

V ^ ^TTTT / 

Certainly, in order to treat the decay of the rj', one must include unitarity 
effects via final state interactions. One very obvious approach is simply to 
include the (energy-dependent) width of the rho-meson in the propagator in the 
vector-dominance form Eq. ^ via 

^-' ^^—W7 -7 (^l) 

TOp - Stttv mi - S^TT - impl p(S7r7i-) 

This use of vector width-modified vector dominance already makes an important 
difference from the simple anomaly—tree level—results (especially in the case 
of the T]'), changing the predicted decay widths from the values 35 eV and 3 
KeV quoted above to the much more realistic numbers 

= 62.3 eV, = 67.5 KeV (52) 
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Figure 5: Shown are contact (a) and VMD (b,c) contributions to ry, ry' ^ tt+tt 7 
decay. 


if the parameters 

Fii/F^ = 1.3, Fo/F^ = 1.04, 6 = -20° (53) 

are employed. However, this procedure does not match onto the one-loop chiral 
form in the low energy limit. 

In order to determine a form for the final state interactions which matches 
onto both the one-loop chiral correction and to the vector dominance result in 
the appropriate limits, we postulate an N/D structure 


y(s, Stttt) — H^_7r7r7(0, 0) 1 C -t- C 


1 -j— 


where Di{s) is the Omnes function and is defined in terms of the p-wave tttt 

phase shifts via 

citefnt 2 


Di{s) = exp 



ds'5{s') 

Uml s'(s' -s-ie) 


(55) 
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and a, c are free parameters to be determined. In order to reproduce the coef¬ 
ficient of the i^(s,r 7 r) function, which contains the rho width, we require c = 1. 
On the other hand, matching onto the VMD result at 0{p^) can be achieved 
by the choice a = lj2m'^. Thus in the case of the rj the form is completely 
determined. Since the rj' spectrum is closely related and is dominated by the 
presence of the rho we shall postulate an identical form for the r]' case. Using 
these forms we can then calculate the decay widths assuming the theoretical 
values for the anomaly. Using the parameters given in Eq. ^ one finds, for 
example, 

ODf P(s) = 65.7 eV, = 66.2KeV 

ii)Dl''^\s) = 69.7 eV, = 77.8KeV (56) 

There is a tendency then for the numbers obtained via the analytic form of the 
Omnes function to be somewhat too high. 



(a) (b) 

Figure 6: Shown is the photon spectrum inrj —> Tr'^Tr~'Y from Gormley et al. |^l| 
as well as various theoretical fits. In the Figure 6a, the dashed line represents 
the (width-modified) VMD model. The (hardly visible) dotted line and the solid 
line represent the final state interaction ansatz Eq.^ with use of the analytic 
and experimental version of the Omnes function respectively. Figure 6b shows 
the experimental Omnes function result (solid line) compared with the one-loop 
result (dotdash line). 

We can also compare the predicted spectra with the corresponding experi¬ 
mentally determined values. As shown in Figure 6, we observe that the experi¬ 
mental spectra are well fit in the rj case in terms of both the N/D or the VMD 
forms, but that the one-loop chiral expression does not provide an adequate 
representation of the data. In the case of the corresponding rj' decay the results 
are shown in Figure 7, wherein we observe that either the unitarized VMD or 
the use of provides a reasonable fit to the data (we get x^/dof=32/I7 

and 20/17, respectively), while the use of the analytic form for the Omnes 
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Figure 7: Shown is the photon spectrum in rj' 7 r+ 7 r “7 from Abele et al.|Q as 
well as various theoretical fits. As in Figure 6 a, the dashed line represents the 
(width-modified) VMD model. The dotted and solid lines represent the final 
state interaction ansatz Eq.^ with use of the analytic and experimental version 
of the Omnes function, respectively. Here the curves have been normalized to 
the same number of events. 


function yields a predicted spectrum (y^/dof=104/17) which is slightly too low 
on the high energy end. However, for both rj and rj' we see that our simple 
ansatz— Eq|54| —provides a very satisfactory representation of the decay spec¬ 
trum. Our conclusion in the previous section was that if the mixing angle and 
pseudoscalar coupling constants were assigned values consistent with present 
theoretical and experimental leanings, then the predicted widths and spectra of 
both 77, ?7 —!■ TT+TT^y are basically consistent with experimental values. Our goal 
in this section is to go the other way, however. That is, using the assumed N/D 
forms for the decay amplitude, and treating the pseudoscalar decay constants 
Eg, Fq as well as the r] — r]' mixing angle 0 as free parameters, we wish to inquire 
as to how well they can be constrained purely from the experimental data on 
?7,77' —> 77 and 77,77' ^ 7r+7r“7 decays, with reasonable assumptions made about 
the final state interaction effects in these two channels. 

On theoretical grounds, one is somewhat more confident about the extraction 
of the threshold amplitude in the case of the lower energy 77 ^ Tr+Tr^y system. 
Indeed, in this case the physical region extends only slightly into the tail of 
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Fg/F, 

Fo/F, 

9 

VMD 

N/D^nal 

N/DfP 

1.28 ±0.24 
1.49 ±0.29 

1.37 ±0.26 

1.07 ±0.48 
1.02 ±0.42 
1.02 ±0.45 

-20.3° ±9.0° 
-22.6° ±9.6° 
-21.2° ±9.3° 


Table 3: Values of the renormalized pseudoscalar coupling constants and the 
r] — rj' mixing angle using the r],!]' — 77 and 77 — tttt^ amplitudes in a three 
parameter fit. 



Fg/F, 

Fo/F, 

9 

VMD 

N/D^nal 

N/DfP 

1.28 ± 0.20 
1.48 ±0.24 
1.38 ±0.22 

1.07 ±0.04 
1.09 ±0.03 
1.06 ±0.03 

-20.8° ±3.2° 
-24.0° ±3.0° 
-22.0° ±3.3° 


Table 4: Values of the renormalized pseudoscalar coupling constants and of the 
77 — 77 ' mixing angle obtained from a maximum likelihood analysis using the 
77 , 77 ' — 77 and 77 , 77 ' — 7 r 7 r 7 amplitudes. 


the rho unlike the related 77' decay wherein the spectrum extends completely 
over the resonance so that there exists considerable sensitivity to details of the 
shape. Thus a first approach might be to utilize only the two-photon decays 
together with the 77 —> 7r+7r“7 width in order to determine the three desired 
parameters. In this fashion one finds the results shown in Table 1. We observe 
that the results are in agreement, both with each other and with the chiral 
symmetry expectations—~ 1.3, Fq/Ftt ~ 1, and 9 ~ —20°. However, 
the uncertainties obtained in this way are uncomfortably high. 

In order to ameliorate this problem, we have also done a maximum likelihood 
fit including the 77' — TTTry decay rate, yielding the results shown in Table 2. We 
observe that the central values stay fixed but that the error bars are somewhat 
reduced. The conclusions are the same, however—substantial renormalization 
for Fg ~ l.dFjr, almost none for Fq ~ F,r, and a mixing angle 9 ~ —20°. These 
numbers appear nearly invariant, regardless of the approach. 

An interesting aside here is the recent observation by Bar and Wiese that 
the TT^ —> 77 reaction alone does not verify the three color hypothesis—in a 
careful analysis the W-dependence of the Wess-Zumino-Witten term is com¬ 
pletely canceled by the W-dependent part of a Goldstone-Wilczek term, and 
that it is only the 77 —> 7r+7r“7 measurement which truly confirms the result 
that Nc = 3@. 

Having above confirmed the basic correctness of the predictions of the 
anomaly (and thereby of this important cornerstone of QCD) we move now 
to the important three pion decay of the eta, which occurs independent of the 
anomaly and which rather probes the conventional two- and four-derivative piece 
of the chiral Lagrangian. 
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2.3 T] —>■ TrTTTr 

The decay of the isoscalar eta to the predominantly 1=1 final state of the three 
pion system occurs primarily on account of the d-u quark mass difference 
and the result arising from lowest order chiral perturbation theory is well- 
known|| 

Amp (778 ^ 7 r“ 7 r^ 7 r“) = S°‘^d‘^^C{sab, Sac, Sbc) + Permutations (57) 


where 

C(s, t, u) 

and we have defined 
Sab = {Pa + Pb)^ 


Bojmd - mu) 
3V3F2 


1 + 


3(s - So) 


m2 — m2 


and 


So = -^{m^ 


■ ml+ + ml- 


mlo). 


(58) 


(59) 


Equivalently, we can write the prefactor of Eq. in a form which respects the 
reparameterization invariance of Kaplan and Manoharp^ 


Boirud - m„) 1 ^ 

- — 

3V3E2 (52 


(60) 


where 

2 rris-m m^ + m 

Q (61) 

rud - m„ md + m„ 

and m = \ {md + m„) is the average u, d quark mass. 

Thus the decay rj ^ Sir can be used in order to determine the quantity Q^. 
Alternatively, if is given from some other process then the 77 —> 37 r amplitude 
is completely determined. The standard approach to such a determination is to 
utilize the pseudoscalar meson masses via the relation 


q 2 ^ rriK rriR " 

m2 (m^o - m^+)QCD 


(62) 


where (m^o — ™^+)qcd is the nonelectromagnetic component of the K^,K~^ 
mass difference— 


(m^o m^+)QCD — (^A'O '^x+)e-xpt ^7i'+)em (^3) 

In order to evaluate the right hand side of Eq. ^ one generally uses Dashen’s 
theorem, which guarantees the identity of the electromagnetic piece of the kaon 
and pion electromagnetic mass shifts in the chiral symmetric limit 

{ml+ - mlo ) = {m\+ - m^o)EM- (64) 


This simple assumption gives then 


Q 


2 _ 
Dashen 


m\ m|. — m 2 

m2 m^o - m^+ + ml^ - mlo 


24.1 


(65) 
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and results in a prediction 


r(77 ^ 7r''"7r 7r°) = 66 eV (66) 

in strong contradiction to the experimental result 

r®"^P(r;^7r+7r-7r°) =281±28 eV. (67) 

At first sight this would appear to be a rather strong and irreparable violation 
of a lowest order chiral prediction and therefore not salvagable by the expected 
0(771^/(471^^)^) ^ 30% corrections from higher order effects. However, this is 
not at all the case. In fact the one-loop and counterterm contributions were 
calculated by Gasser and Leutwyler and were found to enhance the lowest order 
prediction by a factor 2.6, yielding]^ 

ptfteo(^ ^ 7r+7r-7r°) « 167 ± 50 eV, (68) 

which is a significant improvement, but still somewhat too low. The origin of 
such a large correction lies primarily with rjg, 770 mixing which generates a factor 
(cos0 — V^sinfl)^ ~ 2 leaving the expected 30% corrections due to conventional 
higher order loop and counterterm contributions. However, recent work has 
indicated that Eq. is probably a considerable underestimate due a significant 
violation of Dashen’s theorem. Indeed, the Dashen requirement was derived in 
the limit of chiral symmetry —= 0—and plausible estimates of chiral 
breaking effects have yielded the estimate]^ 

Qx-broken ~ O'^QaashenJ Qx-broken ~ 21.7 (69) 

which corresponds to an additional ~40% enhancement of the chiral calculation 
Eq. i.e. r‘^®°(r 7 ^ 7r+7r“7r°) ~ 240 eV, and puts the result in the right ball¬ 
park. There are at least two possible sources for the existence of any remaining 
discrepancy. One is the fact that the estimate for the size of Dashen’s theorem 
violation is just that—an estimate. It is possible that the size of the violation 
is more significant than that given in Eq. leading to an even larger value 
for r(77 ^ 7r+7r“7r°). A second possibility is that the simple one loop estimate 
given in ref. 4 is not sufficient to include the full impact of final state interaction 
effects. This has been demonstrated in other processes where the 1=0 S-wave 
TT — TT plays an important role, as it does hereQ. Indeed the closely related 
K ^ Sir reaction is one such caseQ. 

In order to decide which—if either—possibility obtains it is necessary to 
make careful spectral shape measurements in addition to simple lifetime num¬ 
bers. Also it is necessary to confront such results with precise theoretical calcu¬ 
lations. Phenomenologically, we expand the decay amplitude about the center 
of the Dalitz plot as 

(7(5, t,u) = a[l + Py + 72 /^ -I- -I- ...] (70) 

where 

3(s-so) j VS{t-u) 

y = - 1 — and x= -^—. 71) 

2mrjATj 2mriAri 
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where — 2m^± — is the Q-value. These parameters have been 

determined experimentally to be 


Layter et al.[52] 


Gormley et al.[^ 


Amsler et al. 


154 


P 

P 

P 


0.54 ±0.007 
0.585 ±0.010 
0.470 ±0.075 


7 = 0.017 ±0.014 
7 = 0.105 ±0.015 
7 = 0.055 ±0.135 


(5 = 0.023 ±0.016 
5 = 0.03 ±0.02 
(72) 


to be compared to the one-loop chiral prediction 


P = 0.665 7 = 0.21 S = 0.04 (73) 

Clearly there is general (though certainly not excellent) agreement, suggesting 
the importance of higher order scattering contributions. 

These have been examined by two Swiss collaborations using dispersion re¬ 
lation treatments in order to address the problem of higher order three-body 
scattering effects. The calculation of Anisovich and Leutwyler quotes only 
the integrated decay rate which is is agreement with experiment if the value 
Q = 22.7 ±0.8 is chosen [|^ —consistent with the Dashen theorem violation cal¬ 
culated in . Similarly the integrated decay rate found in the Khuri-Treiman 
calculation by Kambor, Wiesendanger, and Wyler agrees with the experimental 
rate if the value Q = 22.4 ± 0.9 is usedQ. However, these authors also quote 
values for the spectral shape 


P = 0.58 7 = 0.42 (5 = 0.045 
P = 0.58 7 = 0.415 5 = 0.05 (74) 

where the two sets of numbers correspond to different ways of determining the 
experimental input. Obviously the slope parameter P is in general agreement 
with experiment. However, the situation is more complex for the quadratic 
components. In particular the calculated values for 7 , 5 are in good agreement 
with the measurement of Gormley et al. (or Amsler et al.) but not with that 
of Layter et al. However, experimental uncertainties are significant and a new 
high statistics determination of the spectral shape is needed. 

Additional information is available by studying the neutral decay mode 7 —*■ 
37 r°, for which Bose symmetry determines the decay amplitude to be 

Amp (?7 ^ 37 r°) = ± e{x^ ± ?/^)symP (75) 

where 

{x^ ± 2/^)sym = \ XI+ 2/?) = 

i=l 

Here the dispersive calculation by Kambor et al. predicts 

e = -0.028, or e = -0.014 (77) 

depending on how the experimental input is handled. On the experimental 
side the only measurement of the energy dependence until recently had been of 
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limited accuracy 


Amsler et al.[58] : e = —0.044 ±0.046 

Baglin et al.|^ : e = —0.64 ± 0.74 

Abele et al.@ : e = -0.104 ±0.04 (78) 

which is consistent with (but with large experimental uncertainty) the dispersive 
calculation. However, recently a new result of uncprecedented precision has been 
announced from the Crystal Ball group at BNbjs^ 

e =-0.062 ±0.006 ±0.004 (79) 


Clearly this number is significantly larger than expected from the Khuri- 
Treiman calculation, suggesting that new dynamical input is involved. This is 
not unexpected. Indeed the calculation of Kambor et al. utilized the effective 
chiral Lagrangian at 0{p‘^) is input. Clearly from the agreement with experi¬ 
ment this is the main effect, but one also expects contributions from pieces of 
the chiral Lagrangian of 0{p^) such as 


± Ux^)D^^UD^U%{D''UD^U^) 


(80) 


where x = with Bq being a constant and m is the quark mass matrix, 

~ 47rF.n. ~ 1 GeV is the chiral scale. The coefficients of such terms are 
unconstrained by the strictures of chiral invariance and are experimentally un¬ 
determined at present, since they arise at two-loop order. Nevertheless, their 
presence can lead to peices of the decay amplitude of the form 


A ^ Aik ■ qctta ■ qt + A^ik ■ Uaqt ■ qc + k ■ qtqa • Qc) 


A ( Al A2\ 

“±18+T) 

m2 

Pm^^{Ai-A2) 

’2Q^ 8 Q" 

Q y 07 2 (y 

6mr^ 27 


(81) 


If such a dynamical component is present then its size should be set by chiral 
scaling arguments 

iTid-mu 

A2 ±2 

X 

and an isospin relation 


^{dyn) ± J'*' °{dyn) = €°°°{dyn) 


(82) 


must exist between the quadratic parameters for the charged and neutral chan¬ 
nels. Here the symbol dyn indicates the dynamical (i.e., non-rescattering com¬ 
ponent of the coefficient in question and is found by subtracting the theoretical 
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value obtained from the Khuri-Trieman calculation from the experimental quan¬ 
tity. In this way we find, using the Gormley numbers for experimental input, 

7 +-°(dyn) =-0.025 ±0.015 (5+“°(dyn) =-0.02 ± 0.02 (83) 

and 

e°°°{dyn) = -0.034 ± 0.007 or - 0.048 ± 0.007 (84) 

depending on the dynamical input chosen. The comparison 

'y~^~°{dyn) + ^{dyn) = —0.045 ± 0.03 (85) 


vs. 

e°°°((iyn) =-0.034 ±0.007 or -0.045 ±0.007 ( 86 ) 

is obviously satisfactory within errors but again cries out for a high precision 
measurement of the ij —>■ spectrum, such as would be possible using 

WASA. 

Of course, there is one additional test which we can use. Since according to 
isotopic spin invariance the 37 r'^ amplitude at the center of the Dalitz plot must 
be a factor of three larger than the corresponding 7 r+ 7 r“ 7 r° number, the total 
decay rates should differ by the factor 


r(°)(ooo) 32 

r(o)(-H-o) “ 


(87) 


When rescattering corrections are included, the prediction becomes 

r(OOO) _ f 1.43 one loop 

r(-|— 0) \ 1.41 ± 0.03 Khuri-Treiman 


Both numbers are consistent with the value quoted by the Particle Data 
Group @ 


/ r(ooo) 
vr(±-0)j 


1.404 ±0.034 


(89) 


and are in good agreement with the recent Grystal Ball measurement 


/ r(ooo) 

\n+-o)) 


1.44 ±0.09 ±0.01 


(90) 


Glearly there is plenty of challenge in the three-pion sector for an eta facility, 
but there is also interest in examining the remaining radiative modes rj 
7 r*^ 77 , 3717. 


2.4 T] —>■ 7r°77 

For the decay 77 ^ 71^77 chiral symmetry does not play a important role, but 
vector dominance does. This can be seen from the feature that there exists 
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no contribution at all to this process from the tree level two-derivative La- 
grangian. Rather the lowest order chiral contribution arises at one-loop level— 
r] Stt, KKtt 7 r° 77 . However, the Stt intermediate intermediate state is 
suppressed by the factor rud — rUu while the contribution from KKtt is sup¬ 
pressed kinematically. To see this, we define the general decay amplitude 


Amp{r] 7r°77) = D{s, t, u) [ei ■ e 2 qi • 92 - ei ■ <? 2 e 2 • Qi] 
- E{s, t, u) [-61 • e 2 P ■ qip ■ 92 - ei • pe 2 ■ pqi ■ q 2 
+ ei • q 2 e 2 ■ pp ■ qi + ei ■ pe 2 ■ qiP ■ ( 72 ] 

Then from the one-pion-loop contributions from we find 


D^{s, t, u) = ^^^T(s)F(s, ml), E^{s, t,u) =0 

TT 


where T(s) is the lowest order rj —> Stt amplitude given in Eq. ^ and 

with /3(s) = 


sE(s, ml) = 1-1- In 


4ml ,2 f P{s) + 1 


m - 1 


s — Ami 


(91) 


(92) 


(93) 


with a similar expression obtaining for the kaon loop contribution. Calculation 
of the associated rate yields [Q 

Tioop)?? ^ 7’'°77) ~ 4 X 10“^ eV vs. Texp)?? ^ 7r°77) = 0.84 ± 0.18 eV. (94) 


Thus the one-loop chiral contribution plays a very minor role. Consider, how¬ 
ever, the vector dominance diagram shown in Figure 8 , for which 


D{s, t, u) 


2V3 2 ( "^eFlg 


\ m 


V 


p- q 2 - ml p-qi - ml 


my — t 


TOT/ — u 




— cos 6 — 'J 2 -- 

Fs Fq sin 6 


E{s,t,u) = 


2^3 2 f2eFlg 
9 V 


TOT/ — t mir — u 


^ cos 6 — V2^ sin 9 
Fs Fq 


and yielding for the decay rate 

F'vDiv 7 r° 77 ) = 0.31 eV 


(95) 


(96) 


Inclusion of other higher order effects such as the contribution from a pair 
of anomalous terms—TTTTTry and 7771717 —coupled via a pion loop increases this 
estimate to about half the experimental result, but a considerable discrepancy 
remains and should be the focus of future experimental as well as theoretical 
work. 
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A new number from the Brookhaven experiment was first announced at this 
meeting by Nefkens 


^ 7r77) = (0.38 ±0.11) (97) 

is about a factor of two smaller than the Particle Data Group value. The 
problem with the previous measurements is probably associated with eliminating 
the background from other much more probably neutral modes such as 7 —> 
and confirmation using WASA would clearly be welcome. In this regard, spectral 
shape measurements could be helpful, although this will be difficult, since this 
is a low branching ratio—~ 7 x 10“^—experiment. 


2.5 rj —>■ 'KTT'K'y 


The final mode which we shall mention in this report is 7 —> STry, for which on 
the experimental side there exists at present only an upper bound 


r(7? ^ STry) , 
r(77 —>■ TT+TT^TT*^ 


< 0.0024. 


(98) 


Ordinarily the dominant component of a radiative mode such as this is due to 
the the inner bremsstrahlung process, for which the matrix element is 


Amp (77 —>■ 7 r“'" 7 r Tr^y) ~ 


X 


Amp (?7 —> 7 r“'' 7 r 7 r°) 


iee^ 


(2p+ ± k)^ 

{p+ + kf - ml 


(2p_ ± k)^ 

{p- + fc)2 - ml 


(99) 


and it is experimentally difficult to distinguish any direct photon emission. How¬ 
ever, an exception occurs when the non-radiative process is suppressed in some 
fashion, such as occurs, e.g., in the cases of 7r+ —> e+z/gy and 7r+7r“y, 

wherein the nonradiative process is small because of helicity suppression and CP 
violation respectively]^ . One might have anticipated the same enhancement 
mechanism to apply in our case since the nonradiative reaction rj —> 
takes place only due to the relatively small u-d quark mass difference. For exam¬ 
ple, the direct emission associated with the vector-dominance diagrams shown in 
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Figure 9 could be expected to play an important role. However, a careful analy¬ 
sis by D’Ambrosio et al. has shown that this unfortunately does not appear to be 
the casep^. The contribution from this direct—vector-dominated—mechanism 
is found to have the form 

( V 7r+7r-7r°7) = [Pv ’ Po9+- + Pr, ■ P+9-0 + Pr, ■ P-9o+] 

( 100 ) 

where 

Gfj = k ■ p,p'^ - k ■ pjp^ (101) 

and the combination of Bose symmetry and gauge invariance results in a sup¬ 
pression that makes such pieces even smaller than the pion loop component, 
which arise to 0{p‘^) and which are proportional to — irid- In ref. it is 
estimated that the direct emission (DE) component is only a few percent addi¬ 
tion to the inner bremsstrahlung (IB), even at relatively large photon energies 

[(riB-(-DE - riB)/riB ]^^>90 MeV “ X 10“^ (102) 

which appears to be too small to make this a realistic experimental goal. 



Figure 9: Vector dominance diagrams responsible for rj^ —> tt+tt Tr^y. 


3 Conclusions 

We have examined the ”non-rare” decay modes of the eta meson —rj —> 
77 ,7r+7r“7 , Stt, Tr'^yy, Syry —in light of current theoretical knowledge and within 
the general framework of chiral symmetry. While there exist no striking dis¬ 
crepancies observed with respect to any of these predictions, we clearly iden¬ 
tified problem areas wherein additional experimental scrutiny, such as would 
be possible at CELSIUS, would add significantly to our understanding. These 
conclusions can be summarized succinctly as follows: 
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a : 77 —> 77 —a resolution of the discrepancy between the rates measured 
via the Primakoff effect and via QED production is essential to future 
progress in understanding the eta system in general. 

b : ?7 —> 7r+7r“7 —a careful spectral shape measurement would be useful in 
order check the spectral shape with that predicted on fairly solid theoretical 
grounds from the anomalous sector of QCD. 

c : ?7 — > Stt— a precise spectral shape measurement is called for in order 
to determine whether the existing disagreement between experimental findings 
and (what should be) solid theoretical predictions based on chiral perturbation 
theory are cause by an inaccurate value for the d-u quark mass difference or are 
due to the importance of higher order final state interaction effects. 

d : ?7 — > 7r°77—a precision measurement of the Dalitz plot distribution is 
suggested in order to learn the origin of the existing disagreement between the 
experimental rate and that predicted from vector dominance. 

e : ?7 — > Srry—a determination of an actual rate instead of the existing upper 
bound would be of interest, but theoretical indications are that it will be difficult 
to detect other than the inner bremsstrahlung component. 

Clearly there is lots of interesting physics here and a marriage between pre¬ 
cise and solid new experimental data and careful and well-motivated theoretical 
analysis would, I predict, be a happy one, leading to the offspring of a new degree 
of understanding of an important component of low energy phenomenology. 
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